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ABSTRACT 

The classical Blandford & Payne (1982) model for the magnetocentrifugal acceler- 
ation and collimation of a disk-wind is revisited and refined. In the original model, the 
gas is cold and the solution is everywhere subfast magnetosonic. In the present model 
the plasma has a finite temperature and the self-consistent solution of the MHD equa- 
tions starts with a subslow magnetosonic speed which subsequently crosses all critical 
points, at the slow magnetosonic, Alfven and fast magnetosonic separatrix surfaces. 
The superfast magnetosonic solution thus satisfies MHD causality. Downstream of the 
fast magnetosonic critical point the poloidal streamlines overfocus towards the axis 
and the solution is terminated. The validity of the model to disk winds associated 
with young stellar objects is briefly discussed. 

Key words: MHD - plasmas - solar wind - stars: mass loss, atmosphere - ISM: jets 
and outflows - galaxies: jets 



1 INTRODUCTION 

Astrophysical jets are systematically associated with the 
presence of an underlying accretion disk, both observation- 
ally and theoretically (see Konigl & Pudritz 2000 for a recent 
review). In the case of protostellar objects, accretion disks 
are resolved by means of infrared and millimeter surveys 
and interferometric mappings down to scales of a few tens 
of AU. In the optical and the near infrared, HST high res- 
olution images of disks in several jet sources have also been 
obtained (Padgett et al. 1999). With an apparent relation 
found between accretion and ejection in the form of a strong 
correlation between outflow signatures and accretion diag- 
nostics (see e.g. Cabrit et al. 1990, Cabrit & Andre 1991, 
Hartigan et al. 1995), stellar jets seem to be powered by the 
gravitational energy released in the accretion process. 

These facts and considerations have led several authors 
to develop models of disk winds. The pioneering work of 
Bardeen & Berger (1978) on a hydrodynamic radially self- 
similar model of a hot galactic wind was generalized in 
the seminal paper of Blandford & Payne (1982, henceforth 
BP82) by including a rotating magnetic field. In particular, 
in BP82 it was shown that a cold plasma can be launched 
magneto-centrifugally from a Keplerian disk, similarly to a 
bead on a wire, provided that the magnetic field lines are 



sufficiently inclined from the axis. Since then, steady and 
axisymmetric MHD models, self-similar in the radial direc- 
tion, have been successfully analyzed and generalized in the 
literature (see e.g. Contopoulos & Lovelace 1994, henceforth 
CL94, Li 1995, 1996, Ferreira 1997, Ostriker 1997, Vlahakis 
& Tsinganos 1998, henceforth VT98, Lery et al. 1999). 

A major problem is however still open on the validity 
of the various classes of radially self-similar solutions ana- 
lyzed so far. Because, as it is well known since the original 
work of Weber & Davis (1967) on the rotating magnetized 
solar wind in the equatorial region, acceptable outflowing 
solutions must cross smoothly all singularities related to 
the characteristic speeds of the MHD perturbations, i.e., the 
poloidal Alfven velocity and the slow/fast magnetosonic ve- 
locities. However, in radially self-similar equations the crit- 
ical points are not found where the poloidal speed of the 
flow is equal to the characteristic velocities of these mag- 
netosonic waves. In the cold model of BP82 the "modified" 
fast magnetosonic critical point (where t = 1 in the BP82 
notation) is found downstream of the position where the 
poloidal velocity of the wind is equal to the fast magne- 
tosonic velocity. Subsequently it has been shown that this is 
a general property of the axisymmetric steady MHD equa- 
tions: the singularities of the equations coincide with the po- 
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sitions of the limiting characteristics, or separatrices, within 
the hyperbolic domain of the governing equations (Bogov- 
alov 1994, Tsinganos et al. 1996). In particular, Bogovalov 
(1994, 1996) pointed out the key role played by the singu- 
larity occurring at the fast magnetosonic separatrix surface 
(FMSS). Namely, the asymptotic region of the jet is causally 
disconnected from the base of the flow, only for solutions 
that cross the critical point at the FMSS. This means that 
every terminal perturbation or shock does not affect the out- 
flow structure upstream of the position of this critical point. 
And, Tsinganos et al. (1996) have given several analytical 
examples where the true singularities of the equations do 
not coincide with the positions where the governing partial 
differential MHD equations change character from elliptic 
to hyperbolic and vice versa. For the sake of simplicity from 
now on we shall indicate by 'fast/slow magnetosonic singu- 
larity', or in short 'modified fast/slow', the critical points at 
the FMSS/SMSS. 

It turns out that in none of the previous models of disk- 
winds a solution has been found to cross the FMSS. For ex- 
ample, Li (1995, 1996) and Ferreira (1997), starting from the 
accretion disk, succeeded to cross the slow magnetosonic and 
the Alfven ones, but downwind turning points were found 
where the solutions terminate. Such solutions can be con- 
nected to infinity only through a shock, as suggested by 
Gomez de Castro & Pudritz (1993). However in this case, as 
the wind velocity is subfast magnetosonic, a temporal evo- 
lution of the outflow is expected (Ouyed & Pudritz 1997). 

Cylindrically collimated solutions were found by Os- 
triker (1997) for a cold plasma, integrating the MHD system 
upstream from infinity and crossing the Alfven singularity, 
but always in the subfast magnetosonic regime. On the other 
hand, it has been shown that in collimated winds oscillations 
of streamlines are a common feature (Vlahakis & Tsinganos 
1997). It thus seems that cylindrically collimated solutions 
without oscillations correspond to a rather particular choice 
of parameters that completely suppresses such oscillations. 
A slight change in these parameters induces the onset of 
oscillations which increase in amplitude until the configura- 
tion is destroyed (Vlahakis 1998). Since the Ostriker (1997) 
solutions are asymptotically subfast magnetosonic they are 
likely to be sensitive to perturbations from the external 
medium, unlike solutions that really satisfy all the criticality 
conditions. Therefore, such solutions are likely to be struc- 
turally and topologically unstable (Vlahakis 1998). 

However, it has been shown by Contopoulos (1995) that, 
in the restricted case of a purely toroidal magnetic field, a 
smooth crossing of the FMSS is possible. On the other hand 
in such a case an asymptotically cylindrically collimated con- 
figuration is not found; in fact, a new transition to subfast 
magnetosonic velocities must occur anyway for radially self- 
similar winds. The only way out is then to match the su- 
perfast magnetosonic solution with a shock which is in this 
case in the physically disconnected domain. 

In the present study we extend the analysis of BP82, 
CL94 and Contopoulos (1995) showing that an exact and 
simultaneous smooth crossing of all three MHD critical sur- 
faces is possible. In Sec. 2 we define the equations of the hot 
wind in the framework of a radially self-similar approach 
and outline the numerical technique. In Sec. 3 we explore 
the solution topologies in the region around and particularly 
downstream of the FMSS, where the solution terminates, 



while in Sec. 4 are shown the features of a few solutions 
crossing all three critical points with conditions similar to 
those of BP82. Finally, in Sec. 5 we discuss the possible as- 
trophysical applications of these solutions to stellar jets, and 
summarize the main implications of our results in compari- 
son with previous ones obtained by other authors. 



2 MODEL DESCRIPTION 

In order to establish notation, in this Section we give a 
brief derivation of radially self-similar disk- wind models with 
polytropic thermodynamics. The derivation is along the lines 
of a systematic method which unifies all self-similar MHD 
outflows and includes the BP82 model as the simplest case 
(VT98). 



2.1 General definitions and self-similar 
assumption 

In steady (dt = 0) and axisymmetric (d(f> — 0) MHD, the 
poloidal components of the hydromagnetic field (B, V) are 
defined in terms of the magnetic flux function A and mass 
to magnetic flux function ^a(^4) in cylindrical (z ,zu ,(f>) or 
spherical (r ,6 ,<f>) coordinates, as: 



„ A<p 
V x — , 

VJ 



4np 



(1) 



The azimuthal components are defined in terms of the to- 
tal specific angular momentum L(A) and of the corotation 
frequency £l(A), which are functions of A (Tsinganos 1982): 



L(A) =1B (y*-§t.), 0(^ = 1(^-^1* 
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and of the poloidal Alfven number M 



Al 



(3) 

1, V4, and B4, 
(4) 



Trans Alfvenic flows require that, when M 
are finite, i.e.: 

L 2 I A\ _ 1 

— = zu a (A) = zu*a, 

where t*7* is the Alfven cylindrical radius (the Alfven lever 
arm) along the reference field line a = 1, with the dimen- 
sionless variable a defined as some function of the magnetic 
flux function A which can be reversed to give: 



-B* ^1 



A (a) 



(5) 



where B* is a constant with the dimensions of a magnetic 
field. 

As shown in VT98, all existing classes of radially self- 
similar MHD solutions can be constructed by making the 
following two key assumptions: 

(i) the Alfven number M is solely a function of 9, such that 
the Alfven surface is conical: 



M ee M(6) . 



(6) 



(ii) the cylindrical distance vj to the polar axis of some field- 
line labeled by a, normalized to its cylindrical distance vj a 
at the Alfven point is also solely a function of 9: 
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G(0)^. 



(7) 



Following these two assumptions the set of MHD equa- 
tions is reduced to a system of three ordinary differential 
equations in 9 for M(6), G(6) and the 9 -dependence of the 
gas pressure (see VT98 for details). 



2.2 Polytropic thermodynamics 

Depending on the assumptions on the free integrals A(a), 
^a(o), L(a) and ^l(a), a few classes of radially self-similar 
solutions exist (see VT98). For only two of these classes 
a polytropic relationship between the gas pressure and the 
density is admitted: P = Q(a)p 1 , where Q(a) is the spe- 
cific entropy (the first two cases listed in Table 3 of VT98) . 
In such a case A oc a x/2 , *a oc a {x ~ 3/2>/2 , Q oc a" 3/4 , 
L oc a 1 / 4 , and the system of the MHD equations reduces to 
two first order differential equations for M(9) and G(8), sup- 
plemented by the Bernoulli integral which also provides the 
variable tp(8), the angle between a particular poloidal field- 
line and the cylindrical direction zb at the spherical angle 
9. Note that the parameter x (with the same notation as in 
CL94, while in VT98 x was denoted by F) governs the scal- 
ing of the magnetic field, while the rotation law is assumed 
Keplerian. This particular class corresponds to the radially 
self-similar solutions analyzed in CL94 which contains as a 
special case the classical BP82 solution with x = 0.75. 

The full expressions of dM 2 /d0, AG 2 /A9 and ip(0) are 
given in the Appendix, Eqs. (Al) - (A3) The expressions for 
the physical variables become then: 
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V± = _i/4 G 2 -M 2 
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(8) 
(9) 
(10) 

(11) 
(12) 



2.3 Parameters 

At the Alfven radius za+ along the reference field line a = 1, 
we denote by P* and the pressure and density, respec- 
tively. The magnitude of the poloidal magnetic field at this 
Alfven point is — £>* sin 9+fcos + 9+) while the corre- 
sponding poloidal Alfven speed is —V* sin 9+ /cos (i/>* + 6*), 
with B* = 

The expressions of the free integrals defined in Sec. 2.1 
can now be written as: 



A = — — -a ' . 

x 

n = A — a 

E = V 2 ea 



-3/4 



-1/2 



L = \V*m*a /A 
2 _ QM_ 
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(13) 
(14) 

(15) 



where E is the sum of the kinetic, enthalpy, gravitational 
and Poynting energy flux densities per unit of mass flux 
density, 



E(a) 



2 7 - 



QM 



1 P 



r sin 9B& 



(16) 



while Q and A4 are the gravitational constant and the mass 
of the central body, respectively. 

The solution of the system of Eqs. (Al) - (A3) depends 
on the six parameters x, 7, k, A, e and /i, introduced in Eqs. 
( p^| ) - (|l^) (but see the discussion in Sect. 2.4.3 for the free 
parameters of the model) . Note that we have used for the pa- 
rameters a similar but not an identical notation with BP82, 
since it occurred to us that it is better to choose a differ- 
ent normalization. However, in the following we shall outline 
for convenience the correspondence between our parameters 
and those in BP82. 

Let us first discuss the physical meaning of the above 
parameters. First, the exponent x is equal to 3/4 in BP82, 
while in Ferreira (1997) it is related to the ejection index £ = 
2(x — 3/4). This index £ is related to the accretion rate and 
to the mass flux in the wind if also the structure of the disk is 
assumed radially self-similar (see e.g. Ferreira 1997). Second, 
we remind that 7 is the usual polytropic index. Next, the 
constant k is the Keplerian speed at radius vo* on the disk, 
in units of V», i.e., it is proportional to the ratio of the 
Keplerian speed to the poloidal flow speed at the Alfven 
radial distance, V p ,a, and is related to the corresponding 
constant kbp in BP82. Since re is also proportional to the 
mass to magnetic flux ratio, it is often called 'the mass loss 
parameter' (Li 1995, Ferreira 1997), 



™*V 2 



QM 



p,A ■ 



i ty* + 0*) 



= kbp G D 



3/2 



(17) 



The constant A is the specific angular momentum of the flow 
in units of V* ot* and is related to the corresponding constant 
A B p in BP82, 



A = 



L 



V*w* 



AbP^V Go . 



(18) 



The Bernoulli constant e is the sum of the enthalpy, kinetic, 
gravitational and Poynting energy flux densities per unit of 
mass flux density divided by V 2 (along a = 1) and is related 
to the corresponding constant ebp in BP82, 



£bp - 



E 

Finally, the constant p is proportional to the gas entropy. 



(19) 



/'• : 



Mbp [2Go 7-4 sin 



27-2 



27-1 



(20) 



In the above expressions, the label o indicates the respective 
values of G and xj) at the base of the outflow. The correspon- 
dence between the parameter f in BP82 and our ip is 



£0 = cot rj) . 



(21) 



Note that in BP82 7 does not appear since the outflow 
is cold and [i, although it is defined, is never used. Also, a 
similar scaling exists for the parameters used in Li (1995) 
and Ferreira (1997), although with slightly different nota- 
tions and a further relation between x and k (cf. Eq. 28 
in Ferreira 1997) due to the connection with a self-similar 
accretion disk thread by a large scale magnetic field. 
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2.4 Numerical integration 

The numerical solution of Eqs. (A1)-(A3) requires the ful- 
fillment of the regularity conditions at the positions of the 
three singularities (Alfven and slow/fast modified magne- 
tosonic critical points) . This implies that the six parameters 
of the solution are not all independent. In the following we 
first shortly summarize the main properties of the critical 
conditions before we discuss the numerical procedure to ob- 
tain the solutions. 



2.4-1 Critical Points 

It is evident that Eqs. (A2) and (A3) become indeterminate 
at the Alfven surface where G = M = 1. The regularity 
condition at this critical point can be easily found together 
with the value of the derivative of M 2 (see Appendix). Fur- 
thermore, the denominator of Eq. (A2) vanishes when the 
meridional component of the velocity Vg satisfies the quartic 
(Vlahakis 1998): 

Ve - Vi (C 2 3 + Vl) + Clvl iS = , (22) 

where C s is the sound speed, and Va and Va ,e the total and 
meridional components of the Alfven velocity, respectively. 

These singularities are typical 'X-type' critical points, 
and the above equation is the well known dispersion rela- 
tion for the magnetosonic waves. However it is crucial to see 
that these singularities appear not when the flow speed, but 
instead where its meridional component coincides with the 
meridional component of the slow/fast magnetosonic veloc- 
ity. 

Bogovalov (1994, 1996) and Tsinganos et al. (1996) have 
emphasized that the singularities in MHD steady flows do 
not always coincide with the positions where the flow and 
the magnetosonic velocities coincide, but with the limit- 
ing characteristics, i.e., the FMSS and the SMSS. In our 
case the separatrix is found where Ve is equal to either one 
of the triplet of the characteristic speeds (V s ,g, Va,s, V^g). 
This is so because in addition to the azimuthal direction <j> 
due to the assumed axisymmetry, we have a second sym- 
metry direction, which is the radial direction f because of 
the assumed radial self-similarity. Therefore a compressible 
slow/fast MHD wave that preserves those two symmetries 
can only propagate along 8 which is perpendicular to both 
<j> and f; the speed of propagation of such a wave satisfies 
exactly the quartic Eq. ( p2[ ) (for details see Tsinganos et al. 
1996). 

It is obvious that a physically acceptable solution with 
low velocity and high density at the base, but high speed and 
low density asymptotically must smoothly cross at least the 
SMSS and the Alfven singularity. Such solutions have been 
widely analyzed in previous papers and are consistent with 
the observational data on collimated stellar jets. However 
it is unescapable that also the fast magnetosonic singular- 
ity should be regularly crossed in order to have a steady 
structure causally disconnected from the asymptotic region, 
where the jet interacts with the environment (Bogovalov 
1994). 



2-4-2 Numerical technique for the search of solutions 

An inherent difficulty of the problem is due to the fact that 
the positions of the previous critical points are not known 
a-priori, but need to be calculated simultaneously and self- 
consistently with the sought for solution. At these critical 
points we do know some relations between various functions, 
for example, at the Alfven surface the regularity condition, 
Eq. (A4), should be satisfied. However, this knowledge alone 
is not practically enabling us to directly find a solution. 

The way we will follow to construct a solution through 
all critical points is to use the shooting method with suc- 
cessive iterations. By starting the integration from an angle 
8 = 6i we reach a singular point where, e.g., the denomina- 
tor in dM 2 /d8 vanishes, but not the numerator. We then go 
back and change some parameter and integrate again until 
it converges, i.e., the denominator and the numerator van- 
ish simultaneously. A similar procedure is followed to cross 
the other singularity. A rather key point is the choice of 
the starting position of integration. Most of the previous 
studies solved the equations by starting from the equator 
(BP82, CL94), or from infinity (Ostriker 1997). It occurred 
to us that it is more convenient to integrate the equations 
starting from the Alfven critical point, i.e. from the conical 
surface 8 = 9+, and move upstream (towards the base) and 
downstream (towards the external asymptotic region). 

For the numerical integration, besides the parameters, 
we need also to choose the value of the colatitude 8+ and the 
value there of the slope of the square of the Alfven number 
(p* = dM 2 /dO\g t ) together with the angle of expansion of 
the poloidal streamlines (tp*). Some of these quantities must 
be tuned to fulfill the singularity conditions at the three 
critical points. It turned out convenient for the assumed nu- 
merical technique to tune the values of A and p* for getting 
the critical solution. 

Hence, we first prescribe the parameters 7, x, A and k, 
as well as p* , 8+ and ip* while e is deduced from the Bernoulli 
equation, Eq. (A3), and fi from the regularity condition at 
the Alfven point, Eq. (A4). The integration can now start 
upstream from 9 = 8+ and the SMSS is encountered, but we 
cannot pass through it as, e.g., the denominator of dM 2 /d8 
vanishes there, but not the numerator. We integrate again 
with different values of p* until we find the opposite be- 
haviour around the slow magnetosonic singularity (the nu- 
merator vanishes but not the denominator). Iteratively, by 
fine tuning the value of p*, a solution is finally found which 
pass through the SMSS. 

Then we integrate downstream of the Alfven surface 
and the FMSS is encountered, but in general it is not 
crossed. Changing the value of the parameter A we integrate 
upstream again tuning to a new value of p* until the SMSS 
is crossed. Then we integrate downstream towards the fast 
magnetosonic singularity, and repeat all the procedure until 
we find the right values of p* and A that allow the crossing 
of the two singularities. At this point the complete solution 
is obtained by integrating, with the correct values for all the 
parameters, upstream to the base and downstream towards 
the asymptotic region. 
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Figure 1. Solution topologies around the FMSS for 7 = 1.24 (left) and 7 = 1.23 (right). The other parameters are shown in Tab. 1. 
The critical solution drawn with a solid line crosses the line M mj f = 1 at 9 = 8 m j 6° while the noncritical topologies are plotted with 
dashed and dot-dashed lines. For 9 > 8 m j and M m f < 1 the three different branches of the solutions overlap. 



2.4-3 Selection of the parameters and boundary conditions 

In this study, the critical solution depends on the two 'model' 
parameters (7, x) and the three independent 'fieldline' pa- 
rameters (k, 8+, if)*). The remaining ones (e, /i, p+, A) are 
deduced from the Bernoulli equation and the crossing of 
the Alfven, slow magnetosonic and fast magnetosonic singu- 
larities, respectively. This is consistent with the analysis of 
Bogovalov (1997) where it is argued that since the number 
of equations must be equal to the number of independent 
boundary conditions, a unique solution can be found if this 
number of independent boundary conditions equals to the 
number of outgoing waves generated at the reflection of a 
plane wave from that boundary. In f-dependent polytropic 
MHD there are 7 equations and 7 unknowns: the density, the 
pressure, the 3 components of the velocity and the 2 com- 
ponents of the magnetic field. There are also 7 waves: the 
entropy wave and the outwards/inwards propagating slow, 
Alfven and fast MHD waves. So, we need 7 parameters with 
both counts, as expected. Now, if the boundary of the out- 
flow is in the subslow region the number of outgoing waves 
from this boundary is 4, i.e., the entropy, slow, Alfven and 
fast MHD outgoing waves. Subtracting the number of the 
boundary conditions we are left with 3 independent param- 
eters, precisely k, and -1/)*. Note that the polytropic index 
7 and x should not be included in this count since they 
are model parameters and do not depend on a particular 
streamline. 

The integration is terminated in the upstream region 
when M — > M a < 1 and G -> G a < 1, -> O « tt/2 and 
ip = ip . In all the calculations presented here we were able 



to follow the solution up to the equator, i.e. 8 = 7r/2. This 
base should be in principle the disk surface, where our solu- 
tions should consistently fit particular boundary conditions. 
Such an approach has been followed by Ferreira (1997) who 
looked for inflow/outflow MHD solutions with a consistent 
matching on the disk surface (see also Li 1995, 1996). This 
implies some further constraints on the parameters. For in- 
stance, if the disk is also self-similar with a large scale mag- 
netic field, a relation between k and x is expected, i.e., the 
mass loading in the outflow and the magnetic flux distribu- 
tion on the disk. In addition, if the outflow carries away all 
the angular momentum from the disk x must be related to 
A. 

As described above, the procedure to obtain a critical 
solution is extremely lengthy and rather time consuming. We 
must in fact approach as close as possible the singularities 
(A6> w 1(T 3 ), and this requires the determination of the 
parameters up to several digits. As we are mainly concerned 
to analyze the general behaviour of superfast magnetosonic 
solutions, in the present study we do not investigate the 
details of the boundaries of the outflow. Therefore we assume 
that between the base of the wind and the disk surface there 
is a thin 'transition' region that allows the connection of the 
wind with the disk. 

For similar reasons the present analysis has been per- 
formed only for a limited set of values of the parameters. We 
have fixed the 'fieldline' parameters 9* = 59°, ip+ = 40° and 
k = 2, while two values have been assumed for the 'model' 
parameter x: 0.75 (model I) and 0.7525 (model II). In this 
two cases we will assume that the polytropic index 7 < 5/3, 
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i.e. some amount of heating occurs in the plasma. To make a 
comparison with a purely magnetocentrifugally driven out- 
flow, we shortly discuss also the very general properties of an 
adiabatic solution (7 = 5/3) assuming 9+ = 60°, V* = 45°, 
k = 3.873 and x = 0.75 (i.e., values of and i/>* very close to 
those used for the nonadiabatic models I and II on purpose 
have been selected). 

In the following two Sections we outline the main prop- 
erties of the topologies around the FMSS and discuss the 
structure of the critical solutions. 



3 SOLUTION TOPOLOGIES 

We present here the topology of two solutions around the 
fast magnetosonic point, assuming 7 = 1.24 and 7 = 1.23 
for fixed x = 0.75. The two slightly different values of the 
polytropic index define the transition between two families 
of topologies. This drastic change in the topological behav- 
ior of the solutions in the neighborhood of the X-type point 
illustrates the difficulty of exact crossing the fast magne- 
tosonic point. The parameters for the various cases are listed 
in Tab. 1, while in Fig. 1 we plot the two sets of topologies 
for the superfast magnetosonic number M m> f(#). Note that 
this plot is obtained from a projection of the solutions from 
the 3-D space of M(0), G{9) and 9 to the plane M m , f - 9. 
This three-dimensional structure of the topology explains 
why some of the lines obtained by projection are crossing 
each other (see for another such example Tsinganos & Sauty 
1992). This feature of course does not appear in more clas- 
sical topologies of one-dimensional solutions, e.g., Weber & 
Davis (1967). 

In the first case (7 = 1.24) three solutions are plotted in 
Fig. la for different values of A and fj,. The critical solution 
(solid line in Fig. la), moving downstream in the direction 
of decreasing 9 crosses the FMSS at 9m,{ ~ 6°, has a max- 
imum at 6 ~ 0.4° and then at w 0.15° crosses back the 
Afm.f = 1 line but with an infinite slope moving towards in- 
creasing 9. Then, this solution continues marching towards 
increasing 6 and remains always subfast magnetosonic, with 
M m: f reaching a maximum at 9 ~ # m ,f. 

By slightly decreasing A the solution crosses the M m> f = 
1 line with infinite slope at 9 > 9m,{ (dashed line in Fig. la). 
Conversely, for a slightly larger value of A the solution (dot- 
dashed line in Fig. la) reaches a maximum at 9 « # mj f re- 
maining subfast magnetosonic (i.e. it behaves like a 'breeze' 
solution) and becomes superfast magnetosonic with diverg- 
ing slope at 9 « 0.04°. Then, this solution remains always 
in the region 9 < # m ,f, with a spiraling behaviour, i.e., by 
crossing many times up and down the M m> f = 1 transition 
with infinite slope. Note that the solution shown in Fig. 10 
of Ferreira (1997) probably belongs to this family of non 
critical solutions. 

For 7 = 1.23 (Fig. lb) the topology of the non critical 
solutions remains the same. The critical solution however 
shows a different behaviour remaining always in the region 
9 < # m ,f by spiraling around the M m> f = 1 transition (solid 
line in Fig. lb). 

The topological structure of our solutions implies that 
downstream of the FMSS a focal critical point must be 
present, so that no solution can asymptotically reach 9 — 
with superfast magnetosonic speeds. This ought to be ex- 
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7 


P* 


A 2 


M 


Solution 


1.24 


-12.5522 


72.7220 
72.0000 
73.0000 


6.7825 
6.9069 
6.7347 


critical (solid) 
terminated (dashed) 
spiral (dot-dashed) 


1.23 


-12.6468 


75.8919 
75.0000 
77.0000 


6.6983 
6.8506 
6.5092 


critical (solid) 
terminated (dashed) 
spiral (dot-dashed) 


a 


assuming x 


= 0.75, 0* 


= 59°, i/>* 


= 40° and k = 2. 



pected from the construction of this model where we should 
have lime^o V$/Vf,e = 0, if we have a cylindrically colli- 
mated outflow. In other words, the surface M m ,f = 1 needs 
to be crossed again with a downstream superfast /subfast 
magnetosonic transition (see also Contopoulos 1995). At the 
same time, we should keep in mind that these radially self- 
similar solutions are not valid to model outflows around the 
rotational axis, because of their singular behaviour there. 

Note that not all solutions with M m j > 1 are physi- 
cally acceptable because they become subfast magnetosonic, 
crossing the singularity with diverging slope and therefore 
they are multivalued for the same 9. Hence, these solutions 
could correspond to the terminated solutions in Parker's 
terminology for the solar wind with one (Parker 1958), or, 
multiple critical points (Habbal & Tsinganos 1983). Nev- 
ertheless, the present critical solutions are causally discon- 
nected from the inner region of the flow, so that they could 
be stopped by suitable boundary conditions, e.g. through a 
shock with the external medium at some angle 9 min < m ,f 
without affecting the structure of the outflow upstream of 
the FMSS. 

It is worth to mention that, from the technical point 
of view, the main difficulty in obtaining a critical solution 
is the fact that all solutions (critical ones as well as non 
critical ones) always reach M m ,f = 1 with infinite slope at 
some angle 9. They become "terminated" at this point, even 
if they belong to the family of the dot-dashed solution family 
of Figs. 1. And, both families of non critical solutions almost 
coincide far from the vicinity of the critical X-point. This is 
the reason why the crossing of the critical point is so difficult. 



4 RESULTS 

The values of the parameters in the previous Section were 
chosen such as to illustrate the topology of the solution 
around the fast critical point. However, they do not cor- 
respond to some interesting critical solution from the astro- 
physical point of view. For example, the fast magnetosonic 
transition occurs for a rather slow velocity and not far from 
the Alfven critical surface. We found that much more in- 
teresting results are obtained for a flow closer to isothermal 
conditions. We then discuss in the following the properties 
of solutions obtained with 7 = 1.05 and for two sets of the 
remaining parameters (models I and II in Tab. 2). 

In both cases 9+ = 59°, V* = 40° and k = 2, as in the 
previous topological analysis, with x = 0.75 and x = 0.7525. 
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model 


■ 


P* 


A 2 




e 


I 


0.75 


-14.07 


136.9232 


2.9902 


156.617 


II 


0.7525 


-14.02 


136.2261 


3.1715 


158.233 


b assuming 7 


1.05, 0* 


= 59°, V* = 


40° and 


K = 2. 




Figure 2. Plot of the various Mach numbers vs. the vertical 
height z in units of the equatorial cylindrical radius ro = ro(z = 
0) of a particular fieldline. A polytropic radially self-similar model 
is used with the parameters of model I (Tab. 2). The critical tran- 
sitions at the three singularities are marked with vertical lines. 



The remaining parameters are deduced from the require- 
ment to fulfill the criticality conditions and are listed in 
Tab. 2. We remark that this different choice on the scaling 
of the magnetic field x is important to connect the solution 
to an accretion disk in the spirit of what has been done 
by Li (1996) and Ferreira (1997). In such value 
of x larger than some minimum above the value of BP82, 
x = 0.75, is necessary to allow ejection (Ferreira & Pelletier 
1995). However, it does not mean that our solution fulfills all 
requirements to connect to such a disk, as we discuss later. 
The main goal here is to show that the solution is not af- 
fected qualitatively by the change in x as far as the crossing 
of critical points is concerned. 

In Fig. 2 we plot the various Mach numbers along each 
field line a vs. the vertical height z in units of the equatorial 
cylindrical radius w(z = 0) of a particular fieldline and for 
model I. The various critical transitions are indicated, and 
on the disk surface we find G ~ 0.16 and M w 0.02. The 
SMSS almost coincides with the point where the flow be- 
comes superslow magnetosonic (M s , m ~ M B = 1), at z » 0.5. 
The Alfven critical point (M — 1) is crossed at z ~ 3.5 
while the wind becomes superfast magnetosonic (Mf = 1) at 



z « 20. Much farther away is the FMSS, at z ^ 10 4 . Down- 
stream of this position the various Alfven numbers decrease, 
as expected from the previous topological analysis. 

The turning of the solutions is evident in Fig. 3, where 
the poloidal streamlines together with the characteristics are 
plotted. They cross all critical surfaces, and for 8 < 6 m f the 
fieldlines converge towards the symmetry axis such that the 
conical region with 6 < f5 m ,f is causally disconnected from 
the rest of the domain. The two families of the characteris- 
tics in the hyperbolic domain bounded by the cusp and slow 
surfaces are better seen in Fig. 4 obtained for the adiabatic 
case, with a different set of parameters. One family of char- 
acteristics (black) is tangent to the SMSS at M miS = 1 while 
the other (grey) crosses it. Similarly, in the hyperbolic do- 
main bounded by the cone where Mi = 1 one of the family of 
the characteristics (black) is tangent to the FMSS indicated 
by M m .f = 1 while the other (grey) crosses it. We remind 
that the cusp surface (M c = 1) does not coincide with any 
singularity or typical velocity in the flow. 

The components of the outflow speed along a line a = 
const in units of the initial z-component of the flow speed 
at the disk, are plotted in Fig. 5. The units are choosen in 
order to make a direct comparison of this solution with other 
solutions in the literature (e.g. BP82) . Close to the disk level, 



the escape speed is high, V CE 



440, the initial rotational 



speed is lower, V^, = 110 and of the order of the Keplerian 
speed, Vkc-p ~ 31/^0. The azimuthal speed after some 
increase in the region of corotation, approximately up to the 
Alfven critical point at z ~ 4, decays to zero transferring its 
corresponding kinetic energy to poloidal motion. Thus, the 
z- and ^-components of the poloidal motion grow from their 
subslow and subescape values at the disk level where V z = 
1 to the high values obtained at the modified fast critical 
point where V z ~ 10 3 . The poloidal speed exceeds the local 
escape speed around the Alfven transition. A comparison of 
model I and II makes clear that the different values of x 
do not strongly affect the global behaviour of the solutions, 
even though the boundary conditions of the disk are rather 
different. 

Downstream of the Alfven transition the azimuthal 
component of the magnetic field grows to very high values 
in comparison to the poloidal component (Fig. 6). At the 
modified fast critical point practically all the magnetic flux 
is in the azimuthal direction. For example, B^/Bp m 1 at 
the disk, while B^/Bp « 60 after the modified fast tran- 
sition for both, models I and II. From Fig. 6 it may be 
also seen that the flow velocity is largely in the z-direction 
with very small components along ij> and vj. In Fig. 6, the 
main difference between the two models is in the region up- 
stream of the SMSS: for x = 0.75 the angle between the 
poloidal fieldline and the disk surface is ip « 67°, while for 
x — 0.7525 this angle is ip m 56°. Although these values are 
not very different, only the second case matches the outflow 
launching condition for a cold plasma given in BP82. This 
means that magnetocentrifugal driving is more efficient in 
model II at the base. However, we note that at the SSMS 
where the plasma pressure has dropped significantly both 
solutions can be magneto-centrifugally accelerated. The end 
result shown in Fig. 5 is that the terminal speed is lower 
when x is larger, i.e., when the ejection index is higher. This 
result is consistent with Ferreira's (1997) analysis. 

The behaviour of the various components of the con- 
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Figure 3. Poloidal ficldlincs (solid), characteristics (dotted and dot-dashed) and cones of the singular surfaces (dot-dot-dashed) are 
shown for the radially self-similar polytropic solution of model I. In the shadowed regions the governing partial differential equations are 
of elliptic type and no characteristics exist. The parameters are as in Fig. 2. 



served total energy E vs. z, plotted in Fig. 7, provides 
information on the different driving mechanisms that gov- 
ern the dynamics of the outflow. Upstream of the SMSS and 
close to the base, most of the energy flux is electromagnetic 
plus some amount of enthalpy. The kinetic energy of the 
plasma is negligible. As the slow magnetosonic surface is 
approached, the kinetic energy sharply increases with a cor- 
responding decrease of the thermal energy. Downstream of 
the Alfven surface the Poynting flux rapidly decreases; the 
poloidal kinetic energy keeps increasing, becoming largely 
the main component of the energy flux at the position of 
the FMSS. This behaviour is basically the same for both 
models I and II. In order words, there is some contribu- 
tion to the acceleration of thermal origin up to the modified 
slow critical point after which the pressure drops to a rather 



constant value while the magnetic pressure maintains con- 
siderably higher values up to the Alfven transition. 

We conclude this section by pointing out that the two 
solutions we have analyzed here correspond to efficient mag- 
netic rotators in the terminology of Bogovalov & Tsinganos 
(1999), since the ratio of the corotational velocity to the 
poloidal Alfven velocity at the Alfven critical surface (the 
parameter a in their notation) has a value greater than unity 
(«2.13). 



5 DISCUSSION 

Before discussing the main physical implications of our re- 
sults, also in connection with those obtained by other au- 
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equator 

Figure 4. Poloidal fieldlines (solid) and characteristics (gray, black) for a radially self-similar polytropic (adiabatic) model with x = 0.75, 
7 = 5/3, 0* = 60°, V* = 45°, k 2 = 15, A* = 10.9239, A 2 = 2.7935, p+ = -5.5744, e = 9.4487. The cones of the different transitions are 
shown, as well as the regions where equations are elliptical (shaded). 

thors, we show that the present solutions are suitable to 
describe the physical properties of astrophysical outflows. 



5.1 Astrophysical applications 

The modeling of a particular astrophysical outflow requires 
first the calculation of all physical quantities from the non 
dimensional parameters characterizing the particular model. 
We will address here this question of calculating some ob- 
servable quantities of disk- winds associated with protostellar 
objects from the parameters of our model. 

We deduce first the ratios of some characteristic speeds 
at the disk level, keeping in mind that from the numerical 
results we have obtained M ~ 0.01 and G ~ 0.1. We will 
refer in the following mainly to the solutions with x = 0.75. 



First, the ratio of the poloidal Alfven and Keplerian 
speeds at the disk level is: 

The poloidal magnetic field which is essential in the launch- 
ing of the outflow is anchored in the disk and its energy 
density is less than the rotational kinetic energy density of 
the disk. Thus, the field is rather weak to brake the rotation 
of the plasma at the disk and it is carried passively around 
by azimuthal rotation. 

Second, the ratio of the sound and initial speeds at the 
disk level is: 

(g) > .m7x (w r(g)' i p r - M . ,M, 
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Figure 5. Plot vs z of the components of the outflow speed for Figure 7. Plot vs. z of the various components of the conserved 

model I (solid) and model II (dashed), with the parameters given total energy E for model I. 

in Tab. 2. 
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Figure 6. Plot vs. z of the ratios of the magnetic components 
Bipl Bp an d poloidal flow speeds V z /V^ for models I and II (solid 
and dashed lines), with the parameters given in Table 2. 



where V a — V z (z — 0). The initial ejection speed is negligible 
in comparison to the thermal speed at the disk, a situation 
similar to a thermally driven wind. 

Next, the ratio of the sound and Keplerian speeds at 
the disk level is: 



i/ 2 (G \^' 2 



1 



0.314.(25) 



Vk cp J o V0.17 K M^- 1] 

We notice that the Keplerian speed is about 3 times higher 



than the thermal speed at the disk. Thus, thermal effects 
cannot inhibit the rotation of the disk. 

Finally, the ratio of the Keplerian and initial speeds at 
the disk level is: 



10" 



M 



300, 



(26) 



i.e., the initial speed is negligible in comparison to the Ke- 
plerian speed. 

In our case the flow speed at the fast critical point is 
about 10 3 the initial speed V a . In agreement with the obser- 
vations we can reasonably assume the terminal speed of the 
outflow to be ~ 400 km s _1 , such that its velocity at the 
base is V = 0.4 km s _1 . 

In principle, radially self-similar models do not have an 
intrinsic scale length; however from the previous estimate 
of the initial speed one allows to calculate the footpoint of 
the reference fieldline a = 1 on the disk. In units of 10 solar 
radii this cylindrical distance vj is: 



= 0.19 x 



Mo 
0.01 



0T\ 3 1 M 
Go) k 2 M 



Vkms" 1 / 



.(27) 



Hence, for a one solar mass star we get w ~ 12. 5Rq. 

It is also interesting to calculate the mass-loss rate A4 V 
in units of 10~ s Mq yr _1 : 



lO-S-Moyr- 1 



0.0386 



»( 



Mo 
0.01 



2 (B z ,o\ 2 
\WG) 



(28) 



10.Ro 
where 

f(a out ,ai 

and 



Vkms- 1 / 



f{a ut, a in ) . 



-3/4 



a; -3/4 



if 0.75 , 



(29) 
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f(a ut,a in ) = ln— if £ = 0.75. (30) 

By assuming Wi„ = vj , vj out ~ 10vj o and B Zi0 = 8Gwe 
have M w /(W~ 8 Mq yr -1 ) £s 1, with a temperature at the 
base of the flow of: 

T oM = 3xf0V (§j) 4 ^ (^) 2 - 8X10*K,(31) 

r„^ t = r ,,„— ^8x 10 3 K. (32) 

We remind that T is not the temperature of the disk as wc 
have assumed a transition layer between the disk surface and 
the base of the flow (see Sec. 2). This region could be rea- 
sonably related to a corona heated by dissipative processes 
in the plasma (e.g. magnetic reconnection, ohmic heating, 
etc.; see e.g. Konigl & Pudritz 2000). 

As the flow corotates roughly up to the Alfven point 
(Fig. 5) the specific angular momentum carried by the wind 
is J w — A4 w flvj^ while the angular momentum that has to 
be extracted locally at the foot point zu of the fieldline in or- 
der that the disk accretes at a rate M a is J a = (l/2)Q,zUoM a 
(Spruit 1996). If the angular momentum carried by the wind 
is a fraction / of J a while 1 — / is the fraction carried away 
by viscous stresses, then the ratio of the mass fluxes in the 
wind and in the accretion flow is 

^L = £4 £0.015, 

taking into account that vj a = 5.8zu for model I. It follows 
that the rate of the outflowing mass is at most of the order 
of 1% of the rate of the accreted mass; and this is achieved 
when the wind carries all the angular momentum of the 
accreted mass. When the outflow carries a smaller fraction 
of the angular momentum of the disk, the mass loss rate in 
the wind is an even smaller fraction of the mass loss rate in 
the wind. In other words, the mass loss rate in the wind is a 
negligible fraction of the accreted mass, despite that the jet 
may carry most of the angular momentum of the accreted 
mass. Similar results are obtained for the case x = 0.7525. 
Therefore, from the above arguments we may conclude that 
from our solutions we deduce for the physical parameters 
values in reasonable agreement with those observed in this 
class of objects. 

Our solution terminates at z/zu w 2x 10 4 , i.e., at « 400 
Astronomical Units (AU) from the central star. At this po- 
sition we could argue that there exists a shock matching the 
solution with the outermost region of the outflow (Gomez de 
Castro & Pudritz 1993). It is well known that bright knots 
are observed on scales of thousands AU along most protostel- 
lar jets. These configurations are shocks that are interpreted 
as originated either by fluid instabilities on the jet surface or 
by temporal variations in the velocity of the outflow (Bur- 
rows et al. 1996, Ray 1996, 1998, Micono et al. 1998, Konigl 
& Pudritz 2000). It could be reasonable to associate the 
terminal shock of our solutions with the inner knots, found 
at distances down to ~ 100 AU from the star. However we 
cannot ignore that these knots are non steady configura- 
tions and move outwards with velocities ~ 100 -r- 200 km 
s _1 (Ray 1996). We could assume that the shock is well up- 
stream of the optical knots: polarimetric radio data on the 
T Tauri object are consistent with the presence of a shock at 
^20-^40 AU from the star (Ray et al. 1997). Alternatively 



the terminal shock could indeed be located approximately 
at the positions of the inner knots, but there the flow looses 
both self- similarity and steadiness. However as we are in the 
superfast magnetosonic regime, the upwind configuration is 
not affected. Only a much more detailed parametric study 
will be able to test these two possibilities. 

5.2 Physical properties of the critical solutions 

The solutions of this model, in particular Fig. 7, illustrate 
nicely the physical process of transferring electromagnetic 
Poynting energy flux and enthalpy to directed kinetic en- 
ergy flux of the flow in order to accelerate a disk wind 
and then form a jet along the symmetry axis of the sys- 
tem. Thus, the analysis of the previous section is interest- 
ing in the sense that it reveals the driving mechanisms of 
the outflow. The poloidal kinetic energy is negligible at the 
disk level. It then increases rather sharply up to the re- 
gion of the SMSS and Alfven surfaces. This increase is at 
the expense of both, the enthalpy and the electromagnetic 
Poynting energy flux (see, Fig. 7). The poloidal velocity is 
directed basically in the radial direction (Figs. 5 and 6), i.e., 
here part of the random thermal energy together with a part 
of the electromagnetic energy are mostly transformed to di- 
rected wind expansion. Downstream of the Alfven surface it 
is mainly the Poynting energy flux that is effectively trans- 
formed into kinetic energy directed along the rotational axis, 
till the FMSS is encountered. After the FMSS, the flow has 
already reached the maximum speed available from the total 
energy E, which is also approximately equal to the initial 
electromagnetic Poynting energy flux. Then, the accelera- 
tion asymptotically stops. Despite the fact that most of the 
acceleration to high speeds is apparently of magnetic origin, 
the role of the polytropic index and thus of the initial ther- 
mal acceleration may not be negligible, in particular in the 
region before the SMSS. For example, in the case of Fig. 4 
where the flow is exactly adiabatic and 7 = 5/3, the critical 
solution achieves only a very small axial component of the 
velocity which is twice the axial velocity on the equatorial 
plane. In the quasi-isothermal case of models I and II where 
7 = 1.05, the maximum velocity is 1000, higher than the 
equatorial one (Fig. 5). As a matter of fact, this last case is 
closer to the one analyzed in Li (1995) and Ferreira (1997) 
where the gas is isothermal up to the first critical surface 
and then it is taken to be cold afterwards, wherein the pres- 
sure has sufficiently dropped. However, another possibility 
is that the low terminal speed obtained in the adiabatic case 
of Fig. 4 could be due to the lower value of the rotation pa- 
rameter A 2 which is ~ 2.8 in the adiabatic case of Fig. 4, 
as opposed to values ~ 137 and 136 in models I and II and 
similarly for the case examined in Ferreira (1997). 

When the gas has reached a high speed along the z-axis, 
its inertia causes it to lag behind the rotation of the field line 
and the field is wound up, as shown in Fig. 6, resulting to a 
highly twisted magnetic field. Consequently, the strong cur- 
vature force of this predominantly azimuthal magnetic field 
towards the z-axis, causes the poloidal field to collimate. Ini- 
tially the field is flaring away from the rotation axis but the 
curvature force bends the poloidal field lines toward the ro- 
tation axis. The azimuthal velocity peaks around the Alfven 
point which is at a height z = 3.5 and a cylindrical distance 
tv — 5.8 times the starting distance vj in model I. Beyond 
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the Alfven point the rotation drops in accordance to angular 
momentum conservation and thus the centrifugal force be- 
comes negligible. Then, the strong inwards curvature force 
of the twisted field, wins, over the weak outwards centrifugal 
force and gas pressure gradient with the result that the lines 
are bent and eventually collapse towards the rotation axis. 

It is interesting that this feature of the collapse of the 
outflow towards the rotation axis which appears in cold mod- 
els (BP82) and models that do not cross the FMSS (Li 1995, 
Ferreira 1997), is also preserved in the present hot model 
where also all critical points are crossed. This result seems to 
indicate the rather dominant role of the magnetic hoop stress 
in radially self-similar models, contrary to what happens in 
meridionally self-similar models wherein the structure be- 
comes asymptotically cylindrical (Trussoni et al. 1997, Sauty 
et al. 1999, Vlahakis & Tsinganos 1999). 

It is worth to clarify for a moment the term "disk- wind" 
that we used in this study. By that term we simply intend to 
indicate that we describe an outflow from a disk-like struc- 
ture accreting onto a central gravitational object. Thus, the 
flow starts at some angle 9 a above or at the equatorial plane 
of the disk, as opposed to a "stellar" wind flow that starts 
radially above or at a spherical or quasi-spherical source. 
Needless to say that a consistent solution of the accreting 
part of the flow would be required for a consistent solution of 
the inflow-outflow structure in the case of a disk-wind. How- 
ever, such a complete undertaking is beyond the scope of the 
present paper which only intends to emphasize the possibil- 
ity to construct complete steady self-similar solutions for the 
wind crossing all critical points. 

To make such a connection between the disk and the 
outflow, in the spirit of BP82, Li (1995) and Ferreira (1997), 
the first step would be to see how our parameters may fall 
into the range of parameters considered by those models. 
For that purpose, in Eqs. ( p^ ) - (^l]) we have made a corre- 
spondance between our parameters and those used by BP82. 
Thus, in the "standard" solution of BP82 the parameters 
are: kbp = 0.03, Abp = 30 and £' a — 1.58 corresponding to a 
launching angle of the jet at the disk %j) sa 32° < 60°. In our 
case, we find kbp ~ 0.13, Abp ~ 14.57, for both, model I and 
model II. We also have Cbp = cot^c = 0.425 (ip a = 67°) 
for model I and £obp — cottpo — 0.675 (ip = 56°) for 
model II, in the BP82 notation. We note that the values of 
kbp, Abp are close in BP82 and the present model. However, 
the value of the launching angle ip is > 60° in our model 
I because of the additional thermal driving of the outflow 
at the disk level, contrary to the cold model of BP82 where 
ipo ~ 32° < 60° . In summary, our models I and II occupy in 
the space of kbp and Abp, roughly the same domain as in 
BP82 (cf. Fig. 2 in BP82). The only difference is in the value 
of the launching angle ifi which can be larger in the present 
hot model, as expected. These values are within the range 
of the allowed parameters in the (kbp, Abp) space also in 
the analysis of Li (1995, cf. Fig. 3) provided that the mag- 
netic diffusivity is of order one. Note also that model II with 
x = 0.7525 corresponds to an ejection index in the notation 
of Ferreira (1997) £ = 2x - 3/2 = 0.005. 

5.3 Summary 

In this paper we have extended the classical work of Bland- 
ford and Payne (1982), mainly by showing via examples for 



the first time that a solution passing through all MHD crit- 
ical points can indeed be constructed. 

As is well known, the FMSS plays the role of the MHD 
signal horizon such that in an outflow crossing this MHD 
horizon all perturbations which the outflow may encounter 
are convected downstream by the superfast outflow and so 
the steady state solution is maintained. In other words, the 
outflow interior to the FMSS is causally disconnected and 
protected against any conditions it may encounter in the 
interstellar or intergalactic medium towards which the jet 
propagates after it is launched by magnetocentrifugal forces 
from the surface of an accretion disk. 

Unlike other analytical models which produce asymp- 
totically cylindrically collimated outflows (Sauty & 
Tsinganos 1994, Trussoni et al. 1997, VT98, Sauty et al. 
1999, Vlahakis & Tsinganos 1999), this class of radially self- 
similar models cannot continue to infinity but it has to be 
stopped downstream of the FMSS and matched via a MHD 
shock to a subfast outflow that mixes with the interstellar 
medium (Gomez de Castro & Pudritz 1993). This shock can 
connect the present solutions to some breeze, subAlfven or 
subslow magnetosonic branch perhaps also preserving the 
self-similarity. 

Thus, the main difference here with previous results pre- 
sented in the literature is that the asymptotic part of the 
present solutions is causally disconnected from the source 
and hence any perturbation downstream of the superfast 
transition cannot affect the whole structure of the steady 
outflow. 

This task of matching the present solutions with a 
downstream shock however remains a challenge for future 
studies, together with a (time-consuming) more extended 
parametric analysis and also a correct matching of the ideal 
MHD outflow solutions with an inflow in a non-ideal accre- 
tion disk (Ferreira 1997). 
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APPENDIX 

The two first order differential equations for G(9), M(9) 
governing the present class of solutions are: 



dG 2 _ 2G 2 cosi/> 
~~dT ~ sin 9 cos (tp + 6) ' 

dM 2 _ sin (tp + 9) ( ft 2 sin ( 



(1) 



dO 



cos (tp + 9) 



G 



n(x-2)M 



4-27 



M 4 , 2 , cos^ sing M 4 sin 



A V ( *- 2 \l-M 2 



sin (ip + 6) G 4 ' cos 2 (tp + 9) 

M 2 G 4 - M 2 



2M 4 r __ ns (1-G \ +a2 



cos tp 



G 2 1-M 2 
2M 2 — l") G 4 — M 4 1 



sin^ sin (V> + 6>) G 2 (1 — M 2 ) 



7j u (l~M 2 )M- Jri -2X'=- 



M (1-G 2 



G 2 \ 1 — M 2 



+ 



2 M^9 , , 



1 



M 2 cos 2 (tp + 9) 



In the above two equations tp(9) is given by the Bernoulli 
integral: 



tp — 77 — 9 =p arctan 



G 



M 4 sin 2 9 



2e - 



7/i 



2k sin 9 2 
G 



(G - M ) 
G 2 (1 - M 2 ) 2 



(7 — 1) M 2 (t _1 ) 
^ 1/ 2 



+ 



, 1-G 2 
'l-M 2 



- 1 



with the upper sign corresponding to the outflow case (V r > 
0). 

On the Alfven conical surface for 9 — > 9* we have 



1 - G 2 



2 cos tp* 



1-M 2 )^ p„ sinS* cos (tp* + 9*) ' 

where p+ is the slope of the square of the Alfven number. 
Then from Eq. (A2) we get the following third degree poly- 
nomial for p*: 



(x - 2) (4 A 2 + pi sin 2 9*) tan 2 (tp* + 9*) + 
2 (x - 2) 



p J * sin 2 9* + IX- p* + 8y 



tan 9* 



tan(V>* + 6*) + 



(x-2) (up* + p 2 * sin 2 9* + 4A 2 \— ) + n 2 p 2 * sin 9* 

\ tan* 3 9* ) 



(4) 



= 0. 
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